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' Motivated by the recent success of optical trapping of alkali Bose condensate, 
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CN ' we have studied the superfluid state of optically trapped alkali fermions, which 

can have Cooper pairs with total spin J > 2. In this paper, we shall discuss 

^ ■ 

t^J- ' the general structure of these large spin Cooper pairs and their close relation 
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with singlet Cooper pairs with non-zero orbital angular momentum. We also 
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present the exact solution for the J = 2 pairing which shows a surprising 



& ' change of ground state as the spin / of the constituent fermion increases. 

The discovery of Bose-Einstein condensation |LJ in atomic gases has stimulated many new 



research directions. Among these is the search of the superfluid phases of alkali fermions. 
This search has become even more exciting in view of the recent success of confining Bose 



^ ■ condensates in optical traps . Since optical traps are non- magnetic, the spin of the trapped 



atoms are no longer frozen as they were in magnetic traps. This leads to a new class of 
superfluid phenomena. In the case of spin-1 Bose gas like 23 Na and 87 Rb, one of us has 
recently pointed out 23 Na and 87 Rb should have a non-magnetic and ferromagnetic spinor 
condensate respectively, according to the current estimates of their scattering lengths [[|. 
Very recently, experiments at MIT || have verified the basically non-magnetic spinor nature 
of 23 Na and found that its magnetic interaction is indeed antiferromagnetic [R|. 

The physics of alkali fermions in optical traps is equally rich. The fact that all alkali 
fermions (except 6 Li) have hyperfme spins (or simply "spins") / > 1/2 in their lowest 
hyperfme manifold implies that their Cooper pairs can have total spins J > 1. Fermions 
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like 22 Na and 134 Cs which have / = 5/2 and 7/2 can have Cooper pairs with total spin as 
high as 4 and 6. From the example of superfluid 3 He, one can be sure that the internal 
structure of these large spin Cooper pairs will generate a multitude of macroscopic quantum 
phenomena. The purpose of this paper is to point out the structure these large spin Cooper 
pairs, and a surprising change in behavior of a spin- J Cooper pair as a function of fermion 
spin /. 

As a first step, we shall focus on homogenous dilute Fermi gases in zero magnetic fields. 
It is important to understand the homogeneous situation before studying the trapped cases 
||. Moreover, the physics of homogeneous systems are important in their own right. At 
first sight, the weak field limit seems difficult to achieve, for even the Earth magnetic field 
amounts to 10 _4 K, enough to polarize the whole gas. Despite this "strong" background 
field, which can be shielded off to a large extent, one can reduce it effectively to the weak 
field limit by specifying the total spin S of the system. Since the dynamics of these systems 
is spin conserving ||, a prepared spin S out of equilibrium with an external field cannot 
relax to its equilibrium value. The system therefore sees an effective field which would have 
been in equilibrium with the prepared S. By choosing S appropriately, the effective field 
can be made much smaller than the external one. This method has very recently been used 
by Ketterle's group to study the spinor nature of the 23 Na condensate ||. 

Of course, for a pairing state to be observable, its pairing interaction has to be sufficiently 
negative to produce an observable T c . While the scattering lengths of some alkali fermions 
have been calculated, they remain unknown for many alkalis. (See later discussion). In view 
of the lack of information, we have performed a general study of the large spin Cooper pairs. 
In particular, we shall discuss the J = 2 pairing in detail. This is the simplest among all 
large spin pairing which also has an exact solution. The phenomena contained in this case 
reveal the rich physics installed, which turns out to be remarkable indeed. For simplicity, we 
shall call the S-wave spin- J Cooper pairs (made up of two spin / fermions) "spin" Cooper 
pairs, and singlet Cooper pairs with orbital angular momentum J (made up of two spin-1/2 
fermions) "orbital" Cooper pairs. Let us first summarize our findings : 
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(A) The structure of "spin" Cooper pairs is analogous to that of "orbital" Cooper pairs with 
the same angular momentum. This allows one to obtain information of the former from the 
latter, for which an exact solution already exists for J = 2. 

(B) The structure of a spin-J Cooper changes as the spin / of the constituent fermions 
increases beyond a critical value. For Cooper pairs with spin J = 2, they are "ferromagnetic" 
(or "axial") if / > 7/2, but non-magnetic (or "real") if / < 5/2. This change of character 
as a function of / is a result of maximizing the phase space for pairing and is independent 
of interaction parameters, as long as they favor J = 2 pairing. 

Free energy : The low energy effective Hamiltonian of a spin-/ dilute Fermi gas with 
s-wave interactions has been derived in ref. ||. It is rotationally invariant in spin space, 
and is of the form H - fj,N= J dx^+(x)^° (3 (x)^ ( g(x) +\ J dx^+(x)^(x)r a/3; ^^(x)^(x), 
^(x)=-^V 2 5 Q/3 - 7 B-F a/3 , 

2/-1 F 

r^=Es? E < ff<*P\ff; Fm >< //; Fm\fffxu > (1) 

F=0 m=-F 

where M is the mass of the fermion, < / /; Fm\f f \w > is the Clebsch-Gordan coefficient 
for forming a total spin F from two spin-/ particles, gp = iirffap/M, and ap is the s-wave 
scattering length of two spin-/ fermions in the scattering channel with total spin F. Because 
of antisymmetry of the fermions, only even F's appear in eq.([|). 

The order parameter of an S-wave superfluids is ^/^(x) =< ■0 a (x)'0 i a(x) >, which is a 
(2/ + 1) x (2/ + 1) antisymmetric matrix in spin space. For homogeneous systems, is 
independent of x. It is convenient to define the gap function 

Aq,/? = V a p.^ u ^ ^u. (2) 

The free energy according BCS theory is 

1 hr>T °° 1 r if 

T = -TrA+r^A - £ £ Tr [AG(k^ n ) A+G(ko; n )] (3) 

where T _1 is given by eq.([I|) with gp replaced by gp 1 , u n = (2n + l^ksT are the Matsub- 
ara frequencies, G a) g(k, oj n ) and G a p(k, oj n ) are normal Greens functions satisfying matrix 
equation (iu n — H (k)) G(ku n ) = 1, (iu n + H,^ (k)j G{k.uj n ) = 1. 
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General structure of Cooper pairs with spin angular momentum J : Under a spin rotation 
U = exp (-iO- F), ^ a ^{Utp) a . This implies * -> = UAU T , and hence A -> A' = 
UAU T . For gap functions that transform like an angular momentum state |Jm >, they 
must satisfy 



It is easy to see that the solution of eq. (f|) is (Aj^) ft( g oc< ffa(3\ff; Jm >. The general 

structure of the spin-J gap function is therefore 

j j 
(A (J) ) Q/3 oc £ c m < ffa(3\ff;Jm>, or |A (J) >oc £ c m | Jm >, (5) 

m=—J m=—J 

where the second expression in eq.([|) is simply the first written in abstract form. 

To find the that minimizes the energy, and to illustrate the relation of "spin" and 
"orbital" Cooper pairs, it is useful to consider a different representation of A^ J ^. First, we 
note that the singlet state (A (0) ) a/3 oc 77 =< ff;00\ffa(3 > ^2f + 1 satisfies Ur]U T = 77, 
and has the properties r] + r] = 1, U + r] = i]U T , and Ftf = —rjF^. Defining A = S77, 
Eq.(|) then becomes UE$U + = D^E^) , which has the solution [s^] oc [Y Jm (F)] a ^ 
where Yj m (F) is a matrix obtained by first writing the spherical harmonic k J Y^(k) in a 
symmetric rectangular form, and then by replacing by the matrix j7|. For example, 
since k 2 Y 2 i(k) oc k z (k x + ik y ), we have Y 2 i(F) oc F Z (F X + iF y ) + (F x + iF y )F z . The general 
form of the order parameter within the angular momentum J subspace is then 

j 



AJ3= E c m [Y Jm (F)vU- (6) 

m=—J 



Using Wigner-Eckart theorem, it is easily seen that the two representation eq. (|5|) and eq. (0) 
are identical. 

Next, we note that r J Yj m (r) is a homogenous polynomial of r satisfying Laplace's equa- 
tion. It can therefore be written as r J Yj m (r) = A il i 2 ,,,i J ri 1 ri 2 ...r ij , where A^^j is sym- 
metric in all its indices and vanishes whenever any two indices contract. We can then write 
A( J ) as 



A« = £ A lll2 ... lJ F ll F l2 ...F ljV . (7) 

h—ij 

It is also useful to compare the spin structure of A^ J ^ in eq.(|7|) with the orbital structure 
of the singlet Cooper pairs of spin- 1/2 fermions. The order parameter of the latter is 
A(k) =< Cf(k)c|_(— k) >, where c**(k) creates a spin +1/2 fermion with momentum k at the 
Fermi surface. For pairing with even orbital angular momentum J, A^ J ^(k) = J2 m c m Yj m (k), 
or 

A (k) = Ai x i 2 i J ki l ki 2 ...ki J . (8) 

h—ij 

Comparing eq.(|7|) and eq.(||), one finds that they are almost identical except that the Fj's 
are non-commuting matrices whereas the fcjS are c-numbers. On the other hand, this means 
that these two structures approach each other as / increases, since the spin operator F 
behaves more like a classical vector. 

The general scheme for determining A J and the J = 2 pairing : At temperature is 
lowered, superfluid condensation first takes place at the (even) angular momentum channel 
J with largest negative coupling constant gj. The free energy eq.([|) to the quartic order in 
A( J ) is 

T = -^TrA^+A^ + ^Tr(A( J )A( J ) + ) 2 (9) 

where a = JV(0)ln(T c /T), T c = LUe^e- 1 ^^ ^ l.Ue F e- n ^ 2kF \ aj \\ N(0) is the density 
of state at the fermi surface per spin, (3 = 7C(3)/(87r 2 T c 2 ), and €f and kp are the Fermi 
energy and momentum. To determine A^ J ), we substitute eq.(^) into eq.([|) and find the 
matrix A that minimizes the energy. In the following, we shall present the exact solution for 
for S-wave J = 2 Cooper pairs formed by spin-/ fermions. The solutions of J > 2 Cooper 
pairs will be studied elsewhere for they require much lengthier calculations than the J = 2 
case, which is already lengthy. Our method, however, applies to all J > 2 pairs. 
^From eq.(^|), one can see that T is of the form 

a 
T = — —Aij A* q Ti (FiFjF p F q ) + ^ Aij A* ke A pq A* t Tr (Fj Fj F k F e F p F q F s F t ), (10) 
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After evaluating the traces (see Appendix), we find 



T = -&Tr AA + + Pi 



+ 2 (TiA*A) 2 + /3 3 Tr (A 2 A*' 



(11) 
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where I n = Ei=-/ m ", and 5 = - -^2]- 

Eq.(|TT|) is identical to the free energy of a general d-wave singlet superffuids The mini- 



mization problem of eq.(|Tl|) was solved by Mermin ||. Only three equilibrium phases are 
possible ||: 

(I) "Axial" state : When 3 > -0 l + \/3 1 \„A<x Y 22 (F)r], 

(II) "Cyclic" state : When > 3 > — 6/3i, A oc (F 2 + e 2 ™/ 3 F 2 + e 4 ™/ 3 F 2 )r/, 

(III) "Real" state : When fa < -A0 1 -2\0 X \, A oc [(^(F) + ( 2 (Y 22 (F) + Y 2 ^ 2 (F))]r], where 
Ci and ( 2 are real. 

The portion of the phase diagram in 0i — /3 3 space relevant for our discussion is shown in 
fig.l. Using eq.( p^) to eq. (|i~4D , we note that (0i, 3 ) is in region III for / = | and |, and in 
region I for / > ~. The superffuid is therefore a "real" state for / = | and |, but changes 
to the "axial" state for / > |. 

This change of pairing behavior can be understood as follows. As mentioned be- 
fore, as / — ► 00, the order parameters in eq.(|7|) and eq.(H) become identical, and that 
the energy eq.(|9|) becomes that the weak coupling d-wave superfluid, which has an op- 
timum order parameter Y 22 (k) |J. (This state has "more pairing" Y20 and Y 2 ±\ in the 
sense that its absolute square only has point nodes whereas both |F 2 o(k)| 2 and |F 2 ,±i(k)| 2 
have line nodes). On the other hand, in the most quantum case / = |, there are 
four degenerate Fermi surfaces, labelled by m z = ±|, ±|. The structure of the axial 



and the real state are given by \A axial >= |2;2 >= 3y (||, \ > -|§, \ >), |A rea £ > 
Ci|2;0 > +C 2 (|2;2>+|2,-2>)= % 



3 ^3 _|^3 3 + ,1 ^1 _|^l 1 
2' 2 ^ l 2 '2^y^V | 2'2^ I 2 ' 2 ^ 
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(|§, | > — 1|, | >) + (|-<f"? -?r > — \~%\ ~Y > ) ' wnere the state vectors with integer en- 
tries such as 1 2; > means |||; J = 2, m = >, those with half integer entries such as 
||> "T > means |/ = |, m = | > |/ = |; m = ^- >. One can see that the only two Fermi 
surfaces (m = | and ~) are involved in the pairing in axial state, whereas all four Fermi 
surfaces are involved in the pairing of the real state. Since the real state maximizes the 
amount of pairing, it is favored in this extreme quantum case. As / increases, the number 
of Fermi surfaces appearing in the axial state (i.e the spin state \J = 2,m = 2 >) quickly 
increases. By the time when / reaches ~, the real state no longer has the advantage of 
involving most Fermi surfaces, and the system switches to the axial state, where the spin 
operator F begins to resemble a classical vector. We have thus established Statements (A) 
and (B). 

Observability : The long lived alkali fermions which have / > \ in their low- 
est hyperfme manifold are 22 Na, 40 K, 86 Rb, 132 Cs, 134 Cs, and 136 Cs, which have (/ = 
5/2,9/2,5/2,3/2,7/2,9/2) and lifetimes (2.5 yrs, 10 9 yrs, 18 days, 6 days, 2 yrs, 13 days) 
respectively [|l(J]. According to the recent calculation of Greene, Burke, and Bohn the 
scattering lengths of 40 K are positive, hence unfavorable for pairing. At present, there 
is no information about the scattering lengths of the Cs fermions. On the other hand, 
a A = -65(+40, -120)ob, a 2 = -130(+40, -70) a B , a = -145(+40, -65)a B for 86 Rb; and 
a 4 = -l08(+27, -40)a B , a 2 = -115(+32, -50)a B , a = -117(+34, -55)a B for 22 Na , where 
a B is the Bohr radius and the numbers in the bracket are error bars. 

To estimate T c . we use the value of ftF and ep at the center of the trap. For an anisotropic 
trap with frequencies u± and uj z in the x|/-plane and along z, it is easy to show that kpa 2 = 



(o^XcTa)' ep = l^- 1 -^^") 2 ' w ^ere a± = yh/Mu±, where R is the radius of the cloud in the 
a;y-plane related to the total number of particles N as ^- = f (2/+1) ) 5 with A = oj z /u>±. For 
an isotropic trap (A = 1) with uj ± /2tt = 2000Hz, the expression T C (J = 2) = 1.14e F e~ 7r/2feF|a21 
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gives T c ( J = 2) ~ 1.9 x W~ 8 K for 22 Na with iV = 4 x 10 6 atoms and T C (J = 2) ~ 2.3 x 10~ 7 K 
for 86 Rb with N = 10 6 atoms. Since the lowest temperature reached in current BEC 
experiments is 10~ 9 K, these transition temperatures of fermions (which can be made higher 
by increasing the trap frequency or the anisotropy A) appear to be feasible. Since go is most 
negative, singlet instead of J = 2 pairing will first occur in zero field. This, however, does 
not mean that all higher spin pairing states are non-observable. The singlet spin states can 
be efficiently suppressed in a magnetic field (obtained by specifying the total spin of the 
system as mentioned in the introduction), thereby revealing all other higher spin pairing 



states ||11|| . For length reasons, the effects of spin constraints will be discussed elsewhere. 

Final Remarks: We have shown that the superfluid phenomena of alkali fermions become 
amazingly rich once the spin degrees of freedom are released. Should the current efforts of 
cooling alkali fermions to degenerate limit be successful, transferring the degenerate gas into 
an optical trap will help one to uncover the superfluid phases discussed here. Since 132 Cs, 
134 Cs, and 136 Cs have / = 3/2, 7/2, and 9/2 in their lowest hyperfine multiplet respectively, 
if their scattering lengths turned out to be negative, our result predicts that like 22 Na and 
86 Rb, the ground state of 132 Cs will a "real" state, whereas 134 Cs, and 136 Cs will be an 
"axial" state. 

Appendix: Evaluation of the quartic term in eg. ftLu±) : Denoting X4 = 
A i jAl e A pq A* t Tr(F i F j F k F e F p FgF s F t ), we note that T 4 = V[TrU} , where V = 
AijA* M A pq A* st Q^Qbl d^d^ d£dd~ t ' U = e- ia F e~ lb F e~ tc F e- ld F = e~^ F ', and the subscript 
"o" means a = b = c = d = 0. Expanding U in powers of 9, it is easy to see that X 4 = 
S«=i %n)\ ^ 2n (-^ 2n )o> where Im = Y^m=-f m2n - Next, we note that the relation between 9 
and a, ..,d is independent of /. For spin 1/2 systems, the quantity Q=Tr[e~* a ' CT / 2 e _lb ' CT / 2 

e -ic-*/2 e -id-*/2y Tl ^ =< e ~W-*/2 > can be written ag ^ = Q _ I = £ n=12) ^Zl^e 2n . 

Inverting this relation, we obtained 9 2 as a power series of £, or Q. From this expression, 
we can calculate {T>9 2n ) for n = 2 to 4 by calculating T>Q P for p = 1,2, 3, 4. The latter can 
be easily calculated because they involve only spin 1/2 quantities. Evaluating X4 this way 
gives eq.(|i~2D to eq.fllp. 
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11 



